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Abstract 

In this work we present a scheme for construction of solutions for evolutionary PDEs of some 
polynomial types q' t — P(q,q' x , . . .), where P is a polynomial in a finite number of variables. 
This scheme is a generalization of the existing technique for solution of completely integrable 
PDEs using Inverse Scattering of the Sturm-Liouville differential equation. The KdV equation 

/ 3 / 1 /// 

It = 11x + Tfcn 

is a special case, corresponding to type 1 evolutionary equations. We present a complete 
solution of type 0, and present a KdV hierarchy corresponding to infinite number of polynomial 

evolutionary equations rather for j3 — — J* q(y,t)dy then for q[x,t) itself, of the form 

& = i n b n (p' x ), 

1 3 

where bo = —-/3'xxx + -^{P'x) 2 corresponds to the KdV equation and 4(6 n +i)i = ~i(pn)xxx + 
<ti(/3'xb n )x ■ Soliton solutions (i.e. involving pure exponents only) are presented for each such 
evolutionary equation, demonstrating a "simplicity" of the solutions construction. 
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1 Introduction 

Solution of the KdV equation plays a very significant role in the study of evolutionary Partial 
Differential equations (PDEs). Recall that the equation of the form 

Qt = f(q,q'x, ■ ■ •) 

is called an evolutionary PDE for the function of two variables q(x, t) and a "nice" function /. 
d d 

Here q' x = — q(x,t) and q[ — —q(x,t). Some times the solution q(x,t) is also called a flow. 

ox ox 
In particular case we obtain the Korteweg-de Vries equation 

i 3 / 1 /// , , 

It = ~2 qqx I^W- W 

named after two mathematicians D. J. Korteweg and G. de Vries |DJK95| . There is enormous 
amount of constructive examples of solutions for the KdV equation: using algebraic-geometry 
data [Kri77] . flows on manifolds JLax75, Lax76 , flows on Grassmanians Sat89, BBOO, Har92 
ITT79) . flows on Grassmanians using loop-groups |GS85| . 

Usually the differential equation ((TJ) is studied with q(x,0) defined for all values of a; £ 1, 
belonging to a particular class of functions. This is the so called the initial value problem of 
the KdV equation. A standard technique RB88 , CG67 to solve this equation is using the 
Inverse Scattering theory IMG51 , Fad63 of the Sturm Liouville (SL) differential equation with 
a spectral parameter A: 

- y"x + qy = AV (2) 

In this technique, one associates to the parameter q(x,0), called potential a scattering data 
(or a spectral function p(A)), which differs from class to class. Evolving in time the scattering 
data using a simple differential equation (yielding p(X,t)), one then translates back to the 
corresponding to this data potential q(x,t). It turns out that q(x,t) satisfies the KdV equation 
P (see introduction of [RB88] and [CG67| for the idea). 

Using the idea of P. Lax of so called Lax pair one can show that defining linear operators 

P = -J^ +q{x ,t), B=^ +bl §- + b2 

acting on the space of Schwartz functions, one finds that the Lax condition, for an arbitrary 
function w in the Schwartz space 

(^P-[P,B))w = 

3 3 

is equivalent to 61 = — — q, 62 = — — q' x and the evolution of the function q(x,t) according to 

the KdV equation l|T]l. This explains the choice of the differential equation, we are dealing 
with from this point of view. On the other hand, the same equation arises when one constructs 
solutions of the KdV equation using evolutionary vessels, which we further discuss. 

Using a recently developed idea of a SL vessel 2J, the author was able [Mela] to implement 
the scattering technique for SL equation in a manner, which unifies many previous classical 
works IMG51, Fad63,. There was later discovered a notion of evolutionary vessel Vkclv in 
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[Melb] . generalizing the SL vessel QJ, solving the KdV equation {TJ. The vessel ^Kdv is just 
a collection of operators and speces, where operators satisfy some differential and algebraic 
equations, among which we could mention the Lyapunov equation (|I0p . and a variation of the 
wave equation (|35[) . 

In this work we extend the original notion of the KdV vessel, generating a technique for 
solution of a series of evolutionary PDEs, defined by 

P' t =i n b n (l3' x ), P=±£q(y,t)dy (3) 

1 3 

where bo = — -gP'xxx + ^(P'x) 2 correspond to the KdV equation and A(b n +i)' x — —i(bn)'x'xx + 
ki(f3' x bn)' x . The simplicity of the construction of a solution for © is demonstrated by presenting 
a formula of a soliton in Section \4. 21 

Let us emphasize few points related to the presented results. The recursively defined differ- 
ential equations (J3) are part of the KdV hierarchy, because such evolutions commute between 
them (Lemma [13). We present in this work a technique to study non commuting evolutions 
(Section 13. 1[) but this requires an additional research. Finally notice that it is easy to construct 
a vessel in N < oo variables so that its P(x,ti,t2, ■ ■ ■ , tjv) will satisfy equations Q for each 
variable. For this to hold one has to require that the operators of the vessel will evolve with 
respect to the variable t n as the operators of the n-th KdV vessel (I35|l . (|37[) . There is a hope 
that it will enable to reprove the Witten conjecture [Wit91] using these new ideas. 

2 KdV Vessels 

In this section we present KdV vessels and their properties. Each vessel has four important 
notions, related to it: the tau function, the transfer function, the moments and the differential 
ring R,. All of them are defined in the next sections along with their basic properties. 

2.1 Definition 

We define a set of matrices, which will be frequently referred to. We notice that the theory 
presented here can be modified for other choice of these matrices. For example, there exists a 
definition realizing Non Linear Schrodinger equation (see [Melb I Conclusions and Remarks]). 

Definition 1. SL vessel parameters for the SL equation ([2) are defined as follows 
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Using these parameters we define a KdV vessel. 

Definition 2. A KdV evolutionary vessel, associated to the SL vessel parameters is a collection 
of operators and spaces 

K Kdv = {A,B{x,t),X(x,t);o 1 ,a2,j,'y*^,-t);'H,C 2 ), (4) 
where the operators A,X(x,t) :rl—trl t B(x,t) : C 2 — > H satisfy the following vessel conditions 
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(a.e.): 



dx 






(5) 


dx 




= Bo~ 2 B*, 


(6) 






= iA 7T B > 
ox 


(7) 






= iABa 2 B* - iBo 2 B* A* + iB^B*, 


(8) 






= 7 + a 2 B*X~ 1 Ba 1 - a 1 B*X~ 1 Ba 2 , 


(9) 




AX + XA* 


= -B<j x B*, 


(10) 




X* 


= X. 


(11) 



The operator A may be unbounded, but is assumed to generate a Co semi-group. 

Details on the properties of the vessel ma y be found in |Melb| . Notice that there 

is a standard technique for constructing of such a vessel. Suppose that we are given a 2 x 2 
function of complex variable A, realized BGR90 in the following form 

S{\) = I - BgX^iXI - Ay-Bom, I - identity 2x2 matrix 

for which the operators Bo : C 2 , Xo, A : H — > H satisfy 

X = Xo, AXo + X A* + BocriBl = 0. 

Then follow these steps: 

1. Define B(x) as the solution of ((5JI with the initial value Bo, 

2. Define B(x,t) as the solution of ([7} with the initial value B(x) just constructed. Notice 
that for a diagonal operator A, which is very common in these types of realizations, the 
equation J7J is a wave equation, 

3. Define X{x) = X + /* B(y)a 2 B* (y)dy, 

4. Solve ((8| with the initial value X(x). 

So, the vessel conditions (O,©, J7J and |[8j are satisfied by the construction. The equations 
(|10[) and (|11|) are so called permanent, namely if they hold at t — 0, x = then they hold for all 
x, t and this is true from the choice of Bo, Xo. It worth noticing that the formula for X(x, t) can 
be explicitly written, since the second mixed derivatives, applied to X(x, i) are equal, which 
can be easily checked: 

Theorem 3. Let B(x,t) be a solution of the partial differential equations ([5]), (0. Then the 
formula for X(x, t) is as follows 

X(x,t) = X + f B(y,t)a 2 B*(x,t)dy+ 
Jo 

+ [ [iAB(0,s)a 2 B*(0,s) - iB(0, s)a 2 B* (0, s)A* + iB(0, s) 7 B* (0, s)]ds. (12) 
Jo 

As a result, it is a simple check that the formula for X(x,t) satisfies {B), J5). The main 
reason why we consider and call the vessel %)Kdv as a KdV vessel is the following 
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Theorem 4. Let VOkiIV be a KdV vessel, then the potential 

B 2 

q(x,t) = -2— lntdetX^XOM)) 

satisfies KdV equation (JTJ: 

/ 3 / 1 /// 

Sec Melb for the proof, but the main idea behind it is that using the equations of the 
vessel tUifdv one can differentiate J* q(y)dy — —2tr(a 2 B*(x,t)X~ 1 (x,t)B(x,t) (see later in the 
text equation ([25])) to verify the KdV equation (0. 

2.2 The transfer function and the tau-function 

By definition, the transfer function of the vessel QJxdv is defined as follows: 

S(X,x,t) = I - B*(x,t)X' 1 (x,t)(\I - Ay 1 B(x,t)a 1 , (13) 

Notice that poles and singularities of 5* with respect to A are determined by A only. Multipli- 
cation by the function S(\, x) maps [M Va] . [Melc] solutions of the input LDE with the spectral 
parameter A 

\<J2u(\, x) — o\— u(X, x) + 7«(A, x) = 0. (14) 
ox 

to solution of the output LDE with the same spectral parameter 

d 

Xa 2 y{\,x) ~ a 1 (x)—y(X,x) + 'y t y(\,x) = 0. (15) 

The function 7, is defined by the Linkage condition ((9J. So, it follows that the function S(A, x) 
realizes a Backlund transformation of the SL equation (see [Melc] for details). As a result the 
transfer function satisfies the following ODE: 

j^S(A, x) = o-\a 2 \ + J*(x))S{\, x) - S(X, x)a-\a 2 \ + 7). (16) 
From the Lyapunov equation (|10|) . using simple calculations one can show that 

S*{-\,x,t)<T l S(\,x,t) = (Ti, (17) 
S{\,x,t)a^ 1 S*(-\,x,t) = cri" 1 - (18) 

Following the ideas presented in [Melll] we define the tau function of the transfer function 
S^A, x) in the following way 

Definition 5. For a given realization 

S(\,x,t) = I - B*(x,t)X^ 1 (x,t)(\I - AY 1 B[x,t)ai 
tau function r(x) is defined as 

T = det(Xo 1 X(x,t)) (19) 
Notice that from formula ()12p it follows that 

Xo 1 X(x,f) = J + X _1 (0) f B(y,t)a 2 B*(x,t)dy+ 
Jo 

+ X0 1 f [iAB(0, s)a 2 B*(0, s) — iB(Q, s)a 2 B*(0, s)A* + iB(0, s) r yB*(0, s)]ds. 
Jo 
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Since 02 has rank 1, this expression is of the form I + T, for a trace class operator T and since 

Xo is an invertible operator, there exists a non trivial interval (of length at least — — j— ) on 

11^0 II 

which X(x) and t(x) are defined. Recall l( '•!>'.)' that a function F(x) from (a, b) into the group 
G (the set of bounded invertible operators on H of the form I + T, for a trace-class operator T) 
is said to be differentiable if F(x) — I is differentiable as a map into the trace-class operators. 
In our case, 

-^-(Xo L X(x)) = Xo 1 -^-X(a;) = Xq 1 B(x)a 2 B*(x) 
ax ax 

exists in trace-class norm. Israel Gohberg and Mark Krein [IG69 formula 1.14 on p. 163] 
proved that if X^" 1 X(a;) is a differentiable function into G, then t(x) = sp(X,y ] 
differentiable map into C* with 



))0 



T — = s P ((Xq 1 X(a;))" 1 -^(Xo 'X(x))) = sp(X(x)'X" 1 (x)) = 

= ap(B(x)cr 2 B*(x)X~ 1 (x)) = tr(<7 2 B* (x)X' 1 (x)B(x)). (20) 

Since any two realizations of a symmetric function are weakly isomorphic, one obtains from 
standard theorems [BGR90] in realization theory of analytic at infinity functions that they will 
have the same tau function, i.e. this notion is independent of the realization we choose for the 
given function S(X,x). The following theorem was established for the finite dimensional case 
in |Melll) . We present now its generalization for a regular case, i.e. when all the operators of 
the vessel are bounded. 

Theorem 6. For a Sturm Liouville regular vessel the following formula for 7* holds 



Proof: If we denote 



7» = 7 + 



-ft = — , then 7* = 
r 



T 
T 

i 



(21) 



then it is necessary to prove that 



Consider now the formula ([20 



— = tT(a 2 B*(x,t)X' 1 (x,t)B(x,t)). (22) 

T 

Notice that the expression B*{x,t)*K~ 1 (x,t)B{x,i) is the first moment in the expansion of 

" Hid 



S{\,x,t) 



E 



From the general formula of the first moment appearing in [MVbl section 5] (in this article 
the moments are different from the moments in this article by multiplication on a^ 1 from the 
right) 



H (x,t) 



-P 
: P'-P 2 



hf 



(23) 



x sp - stands for the trace in the infinite dimensional space. 



G 



it follows that 



: tr{a 2 H (x)) = -jfl 



as desired. 

As a result, from (|2Up and (|24|1 we obtain that 



(24) 
□ 



/? = -tr((T 2 B*(3;)X _1 (a;)B(a;)) (25) 

2.3 The output differential ring i2* and Moment relations 

Consider now the moments H n (x, i) of the transfer function S(X, x, i), defined from the following 
formula 

H n {x,i) 



-ax, 



(26) 



S(X,x,t) = J - B*X _1 (AJ - Ay 1 Bai = I - ^ 
or more precisely 

Definition 7. TVie n-i/i moment H n (x,t) of a regular vessel 03 is: 

H n (x,t) = B*(x,t)X~ 1 (x,t)A n B(x,t). 

Moments will play a crucial role in the future research and we will need some of their basic 
properties. 

Theorem 8. The following relations between the moments of a KdV vessel ^Budv hold 



H a — Ho = 0, 

H n + (-l) n+1 K + ^(-lfHttnHn^i = 0, 

n-1 

H n + (-l) n+1 H* + Y^i-iT'^HiViK-i-i = 0, 

i=0 

(H n )' x = ai 1 a 2 H„ + i — H n+1 a 2 o-i 1 + a^ 1 'y. t H n — Hnja^ 1 . 



(27) 
(28) 

(29) 
(30) 



Proof: 1. the first equality and (|28l) follow from the symmetry condition (|17|) . The equation 
(J29J) follows from (JTSJ. The last |30]) follows from the differential equation |16J). □ 
Plugging in the formula for (H n )' x the expressions appearing in Definition [1] for the vessel 
parameters, and using (|21[) rewritten in terms of f3, we obtain that 



(H n )' x = 





1 



n + l ~ H n + 1 



' 


1 " 










+ 



-iW'-P 2 ) -13 



H n — H n 





1 



(31) 



Remember that our goal is to construct (solutions of) evolutionary equation so that the deriva- 
tive of /3' t is a polynomial of derivatives of /3 with respect to x. So it is necessary to consider 
the following differential ring of functions: 

Definition 9. The differential ring 7t* is generated by f3 and its derivatives with respect to x. 
The role of the ring 7?., in the simplest case (dim'H < oo) was investigated in [Melllj . 
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2.4 KdV equation in matrix form 



Let us consider in more details the case of KdV vessel. Using vessel conditions one can check 
that the following formula holds: 

|//„ = ^[B*X~ 1 A n B] = (B*)' t X- 1 A n B - B'X^Xpr 1 A n B + B*X~ 1 A n B' t = 
= using evolutionary conditions 0, © 

= (B%(-iA*)X~ 1 A n B - B*X~ 1 (iAX' x - iX' x A* + iBjB*)K~ 1 A n B + B'X' 1 A n (iA)B' x = 

= using ©, © and JlO]) 

d d 
i—H n+ i + i — [H ]aiH n , 



dx 



Consequently, 



d d d 

— H„ =i — H n+1 + i — [H ]aiH„ 
at ox ox 



As a result we obtain that the following PDE holds 



^S{\, x, t) = i^S(\, x, t) + ^[Ho^SiX, x, t). 



(32) 



(33) 



Notice that —[Ho] has entries, which depend on /3(x,t) and its x-derivatives only: 



Ox 



Wo] 



d 



-A 



i{P' xx -2PI3' x ) 



(I 



dx 



and we will see later that — h^ 1 is in . 

dx 

The KdV equation ([1]) arises also as the compatibility condition for the equality of mixed 
derivatives, applied to S(X,x,t). More precisely from 

9 .8 _„ ., d , d _., 

Tx \ Wt ^,t)] = Wt [ Tx Si\ X M 

using 1)161) , (1331) . (|26l) and considering the constants (not depending on A coefficients) one 
obtains that the following formula must hold 



(7*)t = -i7*(H )!cCri + iP±{Ho)xx&i + icTi(H )' x 'y*. 



When we consider this equation for the entries of 7* 



-i{[3' x ~l3 2 ) -13 
13 i 



the regular KdV equation (0. For this it is enough to focus on the 1, 2 entry. 



(34) 

this becomes 



3 Different types of evolution 

In the next section we present a general scheme for the evolution of f3(x, t). This idea generates 
flows, which commute and do not commute with the KdV flow. In Section [3] we present an 
infinite number of commuting with KdV flows. 
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3.1 General Evolution 

Suppose that evolution of the operator B(x,t) with respect to t is determined by the following 
formula 

^(Bcti) = A^A'Bm, = A(Bm + ABm x + ... + A n Bm n ), (35) 

i=0 

where rm are 2x2 matrices, which will be determined later. Let us find first a differential 
equation for X so that the permanency of the Lyapunov equation (|10[) holds. Notice that when 
the condition specj4 n spec— A* — is fulfilled the solution X(x, t) of (|10[) is unique and as 
a result the operator X(x, t) is uniquely determined. On the other hand, when the condition 
specA n spec— A* 7^ holds, we have to solve for X using <(©J using initial condition, which 
does satisfy (|10[) , It turns out that an alternative approach may be applied and the formula 
for X will be explicitly obtained. Differentiating Lyapunov equation (|10|l we find the following 
formula 

n n 

AX' t + X' t A* + B' t a x B* + Bai(B*)' t = AX' t + X' t A* + A(J2 A^rmB*) + £ Bm t B* {A*Y A* 

i=0 i=0 
n n 

= A{X' t + J2 A i Bm i B*) + (X' t + J2 BfhiB*{A*y)A* = C 

j=0 i=0 

So, the derivative Xj must be chosen in such a manner that the last equality holds for an 
arbitrary A and B, satisfying (0. Consider special cases for different values of i: 

0: in this case, B' t = ABm and A(X' t + Bm B*) + (X' t + Brh B*) A* =0. A simple solutions 
of this is 

Xj = —BmoB* — — Yo, mo — mo- 
1: Starting from B' t — A(Bmo + ABui\) we obtain that 

A(X' t + BmoB* + ABm-iB*) + (X' t + BmoB* + Bm x A*B*)A* = 
We have to choose, denoting by Yi the new term: 

X' t — —Yo — ABmiB* — Bfh\B* A* = — Yo — Yi, mo = mo, mi = —mi 
2: B' t = A(Bm + ABmi + A 2 Bm 2 ). Then in order to satisfy 

A(Xt+BmoB*+ABm 1 B*+A 2 Bm 2 B*)+(Xt+BmoB*+Bm 1 A*B*+Bm 2 B*(A*) 2 )A* = 0, 
we have to require 

X' t = -Yo - Yi - Y a , 

where Y 2 = A 2 Bm 2 B* - ABm 2 B* A* + Bm 2 B*(A*) 2 for a self-adjoint 2x2 matrix 
m 2 — fh 2 . 

Generalizing these examples, we obtain the following theorem. 

d n 

Theorem 10. Let —(Bai) = A A n Bm n ([T]). Then in order to obtain the permanency of 
ot n= o 

the Lyapunov equation (1101) . it is sufficient to require m n = (— l) n m n , 

71 

Y„ = ^(-l) l A^'Bm„B*(47, (36) 

i=0 

and 



•X)Yi. (37) 
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Proof: Simple calculations show that Y„ = Y„ and 

^Y„ + Y n A* + A n+1 Bm n B* + Brh n B*(A*)™ = 0. 

As a result, 

AX' t + X' t A* + irlBaiB] = 0. 
at 

Thus the permanency condition of the Lyapunov equation holds. □ 
Definition 11. An n-th KdV vessel is the colloection 

VOnKdV = (A,B(x,t),X(x,t);o- 1 ,a 2 ,j,'y*(x,t);H,C' 2 ), 

where the operators act as for the ordinary KdV vessel and satisfy the equations ((5j), (|6}, fSJ, 
(|10p . and the two evolutionary equations (135|) and (|37p . The number n appearing in the 
definition ofVJnKdv ls called the evolutionary type. 

Let us closer investigate the formula for the evolution of a moment H n 

l H » si 8 

= -H-[B*]X- 1 A JV B + B*^-\X- 1 ]A N B + B"X- 1 A N ^-B 
dt l 1 dt l J at 

= a' 1 J2 rhiB* (A*) i+1 X~ 1 A N B + B*X-\J2 Y i )X- 1 A N B + B*X~ 1 A N £ A i+1 BrUitr^ 1 

i=0 i=0 i=0 

= a' 1 £) rh l B*{A*y+ 1 X- 1 A N B + B*X- 1 (J2 Y i )X~ 1 A N B + £ H N+i+1 miCT^ 

i=0 i=0 i=0 

Thus we obtain the following equation: 

8_ 

dt 



H N = erf 1 (A" Y +1 X- 1 A N B + B*X _1 (^] Y,)X~ 1 A N B + ^ HN+i+xnH^ 1 (38) 



£=0 i=0 i=0 

Examining the formulas for Yj and using the Lyapunov equation (|10[) , it is possible to rewrite 
this formula in the following form, after multiplication by o\ on the right: 

— Hno~i = Yj PiHN+icri + ^ HN+i+io-ia^ 1 m i , 

i=0 i-0 

where Pi are non commutative polynomials of the first n moments, their adjoints and the 
constant matrices o\ , 02 , nii , fhi . Consequently one obtains that a formula for the evolution of 
S(X,x,t) is 

~ n 71 

m i=0 i=0 

where Q is a polynomial of degree n in A with coefficients as for Pi, which makes the equality 
to hold for positive powers of A. 
Since 

§ilL S = §t [a ' 1[a2X + 7 * )S ~ 8a ?( a * X + Tfi 

= a^M'tS + a^(a 2 X + 7.)S* - S^r^aaA + 7) 

n n 

= <T,- 1 (7«);5 + ar 1 (<r 2 A + 7.)KE A^iJS - S E X+^rm + Q]- 

n n 

-KE ^Pi)S ~SJ2 A'+V-V, + Q]a^(a 2 X + 7) 

i=0 i=0 
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and 

r) r) r) n n 

mat 8 =irJ& Pi)s - s 5 x ^ mi+Q] 

n n n 

= (E \W*)s + (E XPi)sL -s' x J2 v+VrV, + q' x 

i—0 i=0 i—0 

7i n 

= (E X(Pi)*)S + (E XPJW^faX + J»)S - Sa^(a 2 X + 7)]- 

z=Q i=Q 

n 

-[a^(a 2 X +7.)fiT- Sa^(a 2 \ + 7)] £ A'+VrV, + Qj, 

From here, we will obtain that the equality of mixed derivatives: 

1±S-—-S (39) 
dtdx* ~ dxdt b > 

is equivalent to 

n n 

kr 1 (7.)J+^r 1 (^A+7»)(^A l P ! )]5'+S[^ A 1+ Vr 1 m 1 ]ar 1 (a 2 A+7)+ar 1 ( ( j2A+7»)g-0ar 1 (a 2 A+7) 

n n n 

= [£ A*(Pi)L + (X) A^Otrr^oaA + 7 .)]S + Sa^(a 2 X + 7)^ A^VfV] + (40) 

i=0 i=0 i=0 

In the KdV case notice that Q = and there is a commutative relation 

TL 71 

A I+ Vf 1 m i ]cr{" 1 ((T2A + 7) = cr 1 " 1 (cT2A + 7)[E A* + cr^m,], 

n 

since for the KdV case ^ A l+1 a 1 ~ 1 mi = iAcrf^o^A + 7). Then the condition (|39[) . multiplied 
by the inverse of 5 on the right, is 

n n n 

v^M't + ^ l (nX + 7»)(E A<P = E A *( p )* + (E A'P)^ V2A + 7 »). 

i=0 i=0 i=0 

Here the coefficient of A (i.e. of constant) is the desired evolutionary equation 

v^M't = -ffx lr y*Pb + (Po)'x + Poo-rV* 
coinciding with (|34|) for Po = iai(Ho)' x cri. Notice that we can rewrite the last equality as 

( r 1 - 1 (7«)J = ( J Po)L + [Po,ar 1 7«], (41) 
where [A, B] = AB — BA is the usual notation for the commutator of matrices. In the more 

n 

general case we obtain that the term Q = ^ A'Qi is present and the evolutionary condition is 

i=o 

more complicated. For this we have to consider the coefficients of A at the equation (|40[l : 

°\ (7*)t + <J\ X l*Pn + o- 1 " 1 7,Qo - Qof 1 " 1 7 = (Po)' x + Poo-r 1 7* + (Qo)L 
which is equivalent to 

<rx\l*)'t = (Po + Qo)' x + {Po,*! 1 !*} - ffrS.Oo + Oo^rS. (42) 

The main disadvantage of these formulas is that it is difficult to obtain a polynomial evolu- 
tionary equation (|42[) . Usually the terms Po and Qq involves integrals. We will see in Section 
[4] that it can be overcome for a special choice of the matrices nn. 
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3.2 Evolutionary type zero 

In the first simple case, when B[<t\ — A(Bmo), we can differentiate with respect to t equality 
((25)) to obtain 

-ft = tr[a 2 cr 1 ~ 1 mo-B*A*X" 1 B + mocrf 1 <TaB*Xr 1 AB + a 2 B*X~ 1 BmoB*^ 1 B] = 
= tr[a 2 a- 1 rnoB'- {-X' 1 A - X" 1 BcnB'X'^B + m a~ 1 a 2 B*X- 1 AB+ 

+ tr[a 2 B*li- 1 Bm B*X- 1 B] = 
= tr[(-CT 2 o- 1 _1 mo + moa^ 1 a 2 )B*X~ 1 AB - crafff 1 moB*X _1 BcriB*X~ 1 B+ 

+ tr[a 2 B*X- 1 Bm B*X- 1 B] = 

m i2 — mi 2 —m 22 1 , J mi2 m 2 2 
_ -Hi — tr 
m.22 J L 

+ tr[<T2-Homo-Ho]- 



= tr 



H a!H ] + 



where mo = mo = 



mn mi2 
mi2 m22 



Thus 



= (mi2 - mi 2 ) tr(<T2-H"i) - m 2 2 tr( 

+ trlaa-B'X-^moB'X- 1 ^ 

= (mi2 — mi2) tr(a 2 Hi) — m 22 tr 

// 

. r 

+ ! in 12 I) | //,, 2;- 

T 



1 

-1 



Hi) - [ mi2 m 2 2 ] HqctiHo 

H x ) + /im 22 tr([ J J 1 #o) + 
// 

. r 

+ | m 2 2 1 Ho 1 2t + ti[o- 2 H m H 




) + 



1 
-1 



Since the last line in this expression is in 7?.*, we have to require that 



(mi2 — mi2) tr(a 2 Hi) — m 22 tr( 



1 

-1 



H 1 )+/3m 22 tr( 




1 



Ho) G K*. 



(43) 



Using (|10p it is immediate that the first equality in the next list holds. Then differentiating 
these equalities and using (|31[) we obtain other lines: 



tr(aiHo) = - tr(A + A*) = const G R, 
tx(a 2 Ho) = -P 6 R, 



tr( 



1 

-1 



H )=ip' x -if 6-R*, 



tr(a 2 Hi) -tr( 





t 



Ho)=y8tr( 



1 




#o) --(£1- 4/^/3 + 2/3 3 ) GiJ, 
2i 



tr( 



t(l3'-p 2 )tr( 



ffi)+j9tr( 




1 



Ho) 



(44) 
(45) 

(46) 
(47) 
(48) 



2i/3 -1 
1 

' 1 




Taking into account formulas (|47[) . (|48[) we have to require in (|43[) that mi2 — mi2 = m22 = 
in order to obtain in (143 [) an expression belonging to R„ . Finally, the formula for f3' t becomes 



Ho) + - 6(/3') 2 ) G R* (49) 
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(for real mn, 77112) 
-P' 



tr{a 2 H m H } 



T 
T 



2r 



L 1 \ Ho 
mn raii 

77112 



mn 77712 
77712 

r 





" 1 " 


H 






mn|- ] 

777ll[ — ] 
T 



L 2r J 

t" t' t" t' 
I — 77712 h z — ?77l2 — 
2T t 2r t 



in other words, we obtain the following theorem: 

Theorem 12. Suppose that B' t = ABima^ + mi2ai) for real constants m,mi2- Then the 
evolutionary equation for j3 is: 

P' t = ~mf3 2 . (50) 

Moreover, mo = 777CT2 + miioi is the only possible equation, for which B' t — ABmo and Xj = 
—BmoB* defines an evolutionary equation of the polynomial form. 

Notice that the solution of this equation is 



mt — a(x) 

for an arbitrary function a(x). For J = 0we obtain that /3(s,0) 



which must be equal 



for the beta function of a SL vessel. In other, words, starting with a SL vessel, realizing @(x), 
the corresponding evolutionary vessel of type will have beta equal to 



P(x,t) = 



f3(x) 



mt - 



0(x) 



mt/3(x) 



3.3 Evolutionary KdV type 2 

Suppose that B satisfies the following evolutionary equation 

B' t ai = -A 2 B' x ai = A 3 Ba 2 + A 2 B'y, 

which mean that 7772 = o~%, mi = 7. Then the following theorem holds 
Theorem 13. The following formulas hold 

(Hn)t = —(H n+ 2)' x — (Hi)' x aiH n+ i — {Ho)' x aiHoo-iH n , 
S' t = -\ 2 S' X - [\(Ho)' x o-i + (Hi)' x ai + (Ho)' x aiH ai}S, 



(51) 



(52) 
(53) 



and 7* satisfies the evolutionary equation (|41 f) with Pq = (Hi) x ai + (Ho)' x ctiHo(Ji. 

Remark: Notice that the expression ctiPq is indeed anti-self-adjoint (as 7* is such) because of 
(|28[) for 77 = 1. Differentiating this expression we obtain 



(#1); + (H*)' x + (H )' x aiHo + Hoai(H )' x = 0, 
and from here multiplying by 01 on both sides we will obtain the desired anti-self-adjointness. 
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Proof: The formula or the moments is straightforward, by using the evolutionary condition 
(|5ip and derived from this evolutionary condition JS). From this formula it follows that the 
evolutionary equation of S is 

S' t = —A 2 5^ - [A(ff )!>i + {H 1 )' x o l + (H )' x aiH ai]S + Q(A). 

Comparing the coefficients of all non-negative powers of A, we come to the conclusion that 
Q(X) = and this proves the evolutionary formula for S. Formula for the evolution of 7, is 
immediate from (|41|l , □ 
Thus it remains to investigate the formula, in order to understand what kind of evolution 
is obtained in this manner. We will see later (n = 2 of Theorem I19|) that the evolutionary 
formula obtained in this case is of the form )3' t = P{PxjPxx> ■ ■ •) f° r a polynomial P. 



4 KdV hierarchy 

4.1 Definition and analysis of the Hierarchy 

Let us generalize the construction of Evolutionary KdV type 2. We will suppose that B satisfies 
the following evolutionary equation 

B' t ai = (V=L4)"B>i. (54) 

In other words, since B' x ai = —ABa 2 — By this requirement is equivalent to the definition 

m n = — (V^l) n a 2 , m n -i = — (v^T) n y. 

Notice first that if the operator B corresponds to a multidimensional vessel, i.e. depends on 
N (possibly infinite number of) variables B = B(x,ti,t2, . . . ,iiv), so that the evolution with 
respect to ti is defined by the formula (|54|) with the power rii, then 



(B' H )' t . = (((V=l)Ar*B' x )' t . = ((V=I)A)"*((a/=I)A)"*B. 



xx 



and the same result holds for {B' t .)' t . since ((^T)A) n '((^T)A) n ^ = ((^f-i)A) n ' ((V^TjA)™*. 
Thus we obtain: 

Lemma 14. Evolutions, defined by the formula (|54|l commute (including commutation with 
the derivative with respect to x). 

Proof: we have seen that the commutation holds for the operator B. In a similar man- 
ner, it can be shown for X, for this we have to prove that (using (|35[) . (|36[) and m n — 
-(V=l) n o-2, m n _! = -(A/=l) n 7) 



d_ 

dt 2 



[J2(-iyA ni - i B(V^i) ni o- 2 B*(A*y+ J2 (-iy a" 1 - 1 - 1 B(V^T) ni jB* [A*y] 



d_ 

dti 



n2 — 1 



[^2(-l) i A n2 - i B(y/^I) n3 a 2 B*(A*) i +Y^ {-l) i A n2 ~ 1 ~ l B(V^l) n2 yB*(A*y] = 

i=0 i=0 

Let us calculate the left hand side of this equality. 



d 



^-[T.{-^yA ni ~ 1 B{^l)^a 2 B*{A*y+ Y. (-l)M 

UL2 j-0 i=0 



i \ ni — 1—i 



B(y^T) ni 7-B*(yr) 1 ] 



"1 

£(-l)M" 1 - I [(^TA) n2 Bi( % /^T) ni (T2B* + B{^/^l) ni a 2 (B*)' x {-V^lA*) n2 ]{A*y+ 

i=0 

roi — 1 

+ £ (-l)M" 1 - 1 - I [(^>l)' l =B;(y^T)"i7 J B*+B( % /^T)"i7( J B*)L(-y^TA*)^](^; 

i=0 
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Dividing this expression by (\/— l) ni+ ™ 2 and using formula (|35[) we can obtain after regrouping, 
that this expression is 

E (-lTA ni -'[A n2 B' x a 2 B* + Ba 2 (B*y x (-A*) n2 ](A*y+ 

i=0 

n± — l 

+ E (-lYA^-^lA^B'^B* + B~ i {B*)' x {-A*nA*y = 

i=0 

= E (-l)^" 1 " 4 !-^™ 2 ^!" 1 ^^* + Boio-^B* (-A*y i2 ](A*y+ 

— 1 

+ E {-iyA ni - 1 - l [-A n2+1 Ba 2 a^ 1 jB"' - B^ya' 1 o 2 B* A* {-A*) n2 \(A*y = 

i=0 

= using a 2 <Ji 7 + 7<r{" CT2 = V - 

= -7^TE(-l)M ni - l yl' l2 _Ba 1 B*(yl*) 1 + (-1)™ 1 A n2 Ba^^B* {A*) ni + 
i=0 

+^\Jl{-iyA ni - l Ba 1 B'"{-A'") n ' 2 {A*y - {-l) n2 A ni B 1 a- 1 a 2 B*{A*) n2 . 

n 

Finally, starting from the Lyapunov equation (|10[) and applying to this summation E A n ~ z (A*y 

i=0 

we obtain from 

n 

J2(~iyA n - l [AX + XA* + BaiB*](A*y = 0, 

i=0 

that 

n n rt 

^2(_iy A *-i+iX(A*y +^2(-l) i A"- i X(Ay +1 = — ^2(—l) i A n ~ i B<TiB*(A*) i . 

i=0 i=0 i=0 

Since the first two sums are telescopic, we obtain 

n 

A n+1 X + (-l)"X(4*) n+1 = - ^(-^i^'BdiB*^*) 1 . (55) 

i=0 

Plugging this formula into the last expression, we obtain 



-V^l E {-lyA^-'A^BaxB'iA*) 1 + (-l) ni A n ' 2 Ba 2 cT^ 1 B*(A*) ni + 

i=0 

+^/^Y2(-iyA ni ~ t Ba 1 B*{-A*y i2 (A*y - (-l) n2 A ni B^a^G 2 B*(A*) n2 = 

= V^lA^[A ni+1 X+ {-l) ni X(A*) ni + 1 } + (-l)" 1 A n2 Ba 2 a~ 1 -yB*(A*) ni - 

-^l[A ni+1 X+(-l) ni X(A*) ni+1 ](-A*) n2 - {-l) n2 A ni B 1 a- 1 a 2 B*{A > ) n2 = 
= V^lA n2+ " 1+1 X+V^T(-l) ni A" 2 X(A*) ni+1 + {-l^A^Ba^-^B'iA*)" 1 - 
-V^lA ni+1 X(-A*) n2 - V^l(-l) ni+n2 X{A*) ni+n2 - {-l) n2 A ni B 1 ai 1 a 2 B*{A*) n2 . 

Finally, exchanging here n\ «-> n 2 and subtracting these two formulas we will obtain the 
following expression (notice that the terms depending on m + n 2 are canceled) 



\n 2 _ 



y^T(-l)" 1 A" 2 X(A*) ni + 1 + {-l) ni A n2 Ba 2 a- 1 ^B*{A*) ni - 

— y/~^~lA ni+1 X(—A*) n2 - (-l) n2 A ni B-ya- 1 a 2 B*(A*) 
-|V=T(-l) n2 A™ 1 X(A*)" 2+1 + (-l) n2 A ni Ba 2 a~ 1 jB*(A*) n2 - 

-j=lA n2+1 X{-A*) ni - {-l) ni A n2 Bya~ 1 a 2 B*(A*)" 1 ] = 
= -yf^l{-l) n2 A ni [AX- V^lB'ya-^iB* +XA* - *J=lBo 2 a~ 1 >yB*]{A*) n2 - 
— \Z^l( — l) ni A™ 2 [AX — \T^lB~ftj^ 1 o 2 B* + XA* - ^/^lBa 2 a± 1 -yB*](A*) ni = 
-y^T(-l)" 2 A ni [AX + XA* + Ba 1 B*]{A*) n2 - 
-V^l{-l) ni A" 2 [AX + XA" + BuiB*]^*)" 1 =0-0 = 0. 
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□ 

Under these conditions we will see that the evolutionary equation for S and for f3 will 
correspond to the KdV hierarchy. Let us define a sequence K n of matrix-valued functions as 
follows: 

Definition 15. K n (x,t) is the n-th moments of the following function 

K n (x,t) 



K(X,x,t) = -S' x (\x,t)S~ 1 (\,x,t)a^ 1 = 



i=0 

From this definition we are able to find a recursive formula for K n 's involving the moments 
H n of S(X,x,t). Plugging the formula (|26|) into Definition of K n we obtain 

\n+l ffl v ~ 2—1 \n+l ' ~ 2—1 
i=0 i=0 i=0 



from where we get 



K =H' , (56) 

Ki=H[+H^a 1 H , (57) 

n-l 

K n =H' n + Y J K t a 1 H n - 1 . l . (58) 



Lemma 16. The moments K n satisfy the following equality K n = (y — T) n K„ ■ 
Proof: Plugging S _1 *(-A) from (JTZJ) we find that 

= o-^ 1 (criS(X)a^ 1 )(o-iS(X)a^ 1 )' x 
= -S(X)S- L (X)(S(X))' :t S- 1 (X)a^ 
= -K(X). 

and the result follows. □ 
We will see evantually that the moments K n implement the family of solutions (|54[) by 
choosing Poci = (V~ l) n K n in the formula (|4ip . As a result for the understanding of the 
evolution equation (|41|l . we need to understand the entries of 

K' n — a^ 1 -f,K n + KnjtO-i 1 (to be equal to erf (^Yt^i )■ 

Notice that using S' x appearing in ()16|) 

K — — S x S^ 1 a^ 1 = — CT 1 " 1 ((J2A + 7»)<T 1 f 1 + Sa^ 1 (a2X + -y)S~ 1 ai 1 . 

and plugging this formula for K into K' n — '-f t ,K n +K n "f,a^ 1 we find, denoting T = erf 1 (A<72 + 
7), that 

= — CTj _1 (7*)^(T 1 1 +Aa 1 1 a2SrS' _1 (J 1 1 — XSTS^ 1 cr 1 1 (T20" 1 1 +Aa 1 ' 1 7»a 1 1 <T2<7 1 1 — Xo 1 1 O20-^ 1 ~f t o 1 1 . 

And thus it is necessary to understand the term SYS~ X , since all the other coefficients are in 
7?.». Simple calculations show that 
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Denote the moments of SFS 1 as 



Cln bn 
Cn d n 

find the following relations between these moments 



, and equate the coefficients of powers of A to 



" / 1/ 




O.n + 1 b n +l 




dn+1 ^n+l 




c' d' 


= o 1 1 a- 2 


Cn+1 dn+1 




Cn+1 dn+1 





+ CTl 7» 



or the following recursive relations 



a' n = ~b n +i + ic n + i(f3' ~ {3 2 )b n , 

b' n = 2/3b n + i(d n — a„), 

c' n = (o n +l - dn+l) + ~ /? 2 )(dn - Or, 

d' n = b„+i - ic n - i(fi' - f3 2 )b„. 



2[3c n 



Notice that from the first and the last equations it follows that a n+1 + d'„ +1 — 0, and we will 
require 

Definition 17. (Normalization condition) d n = — a n . 

Notice that this must be checked for the initial value only. Under this condition, the last 
system of equations may be rewritten as follows 

( b n +i =-a' n + ic n +i(f3' ~ /3 2 )b n , 
I 2a n +i = c' n + 2i(f3' - (3 2 )a n + 2(3c„, 
{ b' n = 2f3b n -2ia n . 

We can see here that actually the entries a n +i,6 n +i are recursively defined from the previous 
ones, but the last condition also imposes an additional relation between them: 

Theorem 18. Under normalization condition for 1, 2, . . . , n — 1, the entries of b n 's satisfy the 
following recursive differential relation 



K+l = -»n + 4i(^'6n)'. 



(59) 



Remark: This formula is very similar to the formula of deriving conservative functionals of 
the KdV equation |Lax761 Formula 3.3], [CG74| . 
Proof: Let us rewrite the last equation in the system as 

2/36„+i — 2ia„+i = b n+1 , 

where we plug the recursive formulas for a n +i, b n +\. The we obtain 

2P{-a' n + ic„ + i(/3' - (3 2 )b n ) - i(c' n + 2i(f3' - (3 2 )a„ + 2/3c n ) = ~a' n + ic' n + i[(f3' - f3 2 )b„]' . 

Rearranging the terms and performing some calculations this can be rewritten as 

2ic' n = < - i[(/3' - 2{3 2 )b n ]' - if3b' n + i/3'b' n . 

Plugging here 2a n — i(b' n — 2f3b n ), obtained from the last equation in the system we obtain: 

Aic' n = <' + 4i([f3 2 b n ]' - [3b n - P'b' n - /?"&„). 
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Finally, plugging this last formula for 4ic4 in the derivative of the first formula of the system 
(recursive formula for b n +i), we obtain 

4b' n+1 = -2i(b' n - 2/3b n )" + i6™ + 4i([/3 2 & n ]' - 0&£ - p'b' n - 0"bn) + 4«[(/3' - /3 2 )b n }' = 
= -iC+4i(/3'&„)', 

after cancellations. □ 
It is remained to show that the evolution, defined by (|54[) is implemented by K n . 

Theorem 19. For the evolutionary equation (|54|l of type n the following formulas hold 



(H N )' t = {V^l) n [{H N+n )' x + KocnHN-i + ... K N H N , (60) 
S' t = (V=l\) n S' x + ^T[\ n ~ 1 K l) a 1 + \ n ~ 2 Ki + ...+ K n ]<TiS, (61) 

and as a result, the formula for the evolution of fj is 

<»'r 1 (7*)t =i n {K n y x o-i+i n \Kncn,cri 1 j*], (62) 

1 3 

which can be solved recursively using (|59[) starting from bo — — -/3' xxx + -(/3' x ) 2 . 



Proof: Notice first that the formula (|61[) follows immediately from (I60|l at least for big powers 
of A. It is a simple check that it also holds for the small coefficients too. 
So let us show that the formula (|60p holds. Using (|54p and (|37p 

(H N )' t = (B*X~ 1 A N B)' t = (B')' t X' 1 A N B + B*(X- 1 )' t A N B + B'X- 1 A N B' t = 

= (B*)' x (—>/=lA*) n X- 1 A N B + B*X~ 1 (Y n + Y n - 1 )X~ 1 A N B + B'X' 1 A N (y/^TA) 71 B' x 



Plugging the formula (|36|) for Y„_i, Y n and definitions for m„, m„_i, the last formula becomes 

= (B*)' x (-A*) rl X- 1 A N B- jr{-iyB*X- 1 A n - 1 Ba 2 B*(A*yX- 1 A N B- 
(V 31 )' 1 *=o 

n— 1 

- ^(-l) l B*X- 1 yl n - 1 - I B7B*(yl*) l X- 1 yl JV J B + B*X"M JV J 4 n B;. 
From the formula (|55|) it follows that 

n-l 



i=0 

Plugging this formula into the last expression we obtain: 



("V - i=0 
n— 1 

-^B*X" 1 A(ABa 1 +B 7 )B*(yl*) 8 X"M JV B-(-l) n B*X" 1 Bo-2B*( J 4*) n X~ 1 A JV B+B*X" 1 A JV+n 

j=0 

= (ff JV+ „); + Y t (H n - 1 - i y x (-l) i <7 1 B*(Ayx- 1 A N B, 
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after some regrouping. So, it is remained to prove that 

n — 1 n 

(H N+n )' x + ^2(H n -i-iY x (-iyaiB*(A*yx- l A N B =^K i< 7 1 H N+n - i . 

i=0 i=0 

But this follows from the calculation of the coefficient of powers of A~™ _1 fi~ N ~ 1 at the expres- 
sion 

5i(A)S- 1 (A)S(/i) 

in two ways. First from the definition of K n as the coefficient of S' x (X)S~ 1 (X)cri it follows that 
the coefficient of A _n_1 ^i _JV_1 is 

n 

KiUlHN + n-i- 

i=0 

On the other hand, Since S~ 1 (X)ai = a^ 1 S*(— X)ai and since using (|10[l it is easy to show 
that 

S*(-X)<nS(p) - o-i = (-A + p)aiB*(-A/ - A*)'^' 1 (fxl - A) _1 .Bcti, 
we will obtain the desired result. □ 



4.2 Solitons as examples of solutions 

One of the advantages of this approach is that we are able to create solutions of the n-th KdV 
equation. For example taking a finite dimensional space "H = C 3 and defining the operators as 
follows 









" -ik'i 














A 


= diag(— ik?) — 





-ikl 





= -A*, 


















-ikl 
























" 1 


ik\ 


B(x,t) 







e z 











1 


ik-2 


















1 


ik 3 



{k i +k j )x+(k 2n + 1 + k 2n+1 )t 

X(s,t) =I+[ e - — 6,6*]. 

fW ~r Ivj 



we will obtain that B(x,t) satisfies differential equations ([5]), (|35[) and X(a;,t) satisfies ([6]), (|37[) . 
In other words, t(:e) = det(X(:r,f)) creates a soliton for the n-th KdV equation via formula 
(ED). 

Taking in a different manner, we can obtain other solutions of J3J, similarly to the 
solution of the equation {TJ, presented in Mclb . 
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